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A method for solving few-body scattering equations is proposed and examined. The 
C^l ' solution of the scattering equations at complex energies is analytically continued to get scat- 

tering T-matrix with real positive energy. Numerical examples document that the method 
works well for two-nucleon scattering and three-nucleon scattering, if the set of complex 
energies is properly chosen. 

O . §1. Introduction 

m 
Q 

Scattering by more than two particles upon each other requires boundary condi- 
tions in configuration space of increasing complexity with growing particle numbers. 
Those boundary conditions are reflected in momentum space by e — > limits of differ- 
ent resolvent operators G(E + ie), where E is real and above some thresholds. While 
such a singularity can easily be handled on a two-body fragmentation cut, it causes 
already severe difficulties on a three-body fragmentation cut like in the three-body 
problem. There it leads to so called moving singularities x )' 2 ). They can be treated 
directly for instance using spline interpolations or just straightforward subtraction 
methods 3 ^. They can also be circumvented by contour deformations of the path of 
momentum integration 4 ) ~ 7 ) . In any case this requires great care. For more than 
three particles the four-body break up singularities pose even harder challenges. 

In the past several additional methods have been proposed especially for two- 
body fragmentations. For instance in 8 ) the T- matrix is evaluated at negative ener- 
gies , where bound state methods can be employed, and then analytically continued 
to the energy in the scattering region. In 9 )" n ) the T-matrix is evaluated at several 
complex energies above the positive real energy axis and then the value on the real 
energy axis is again determined by analytical continuation. All the applications we 
are aware of are in the field of atomic physics. 

The Lorentz transformation method used intensively by the Trento group 12 ) 
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solves many body equations based on a complex energy, thereby avoiding the com- 
plicated boundary conditions in configuration space or when used in momentum 
space 13 ) avoiding singularities. But in this case no continuation to the real energy 
axis is performed and the imaginary part of the energy is kept fixed. Certain inver- 
sion procedures are subsequently employed. 

Here we would like to reconsider the approach evaluating the two- and three- 
body T-matrix for complex energies and performing an analytical continuation to 
the real axis. In view of future applications to the field of nuclear physics we shall 
employ nucleon - nucleon forces . Section 2 displays our studies for two nucleons. We 
add a brief example for three-body scattering in section 3 end with a brief outlook 
in section 4. 



§2. Two- Nucleon Scattering 



We start with a simple model for the NN forces allowing for an analytical solu- 
tion, the Yamaguchi separable force 14 -* 



1 1 

y + p 2 p' 2 + p 2 

It leads to the two-body T- matrix 



with 



t(p,p';E)=r(E) * * (2-2) 

p l + p l p' + (3 2 



, „. , „ i 7r(B 2 — k 2 )m nkm , , 

t(E) = (A" 1 + ■ ' + i———)- 1 (2-3) 



Af3{(3 2 + k 2 ) 2 2(f3 2 + k 2 ) 2 - 

Now we evaluate t(E) replacing energy E = k 2 /m by E + ie, where e is now finite 
and we choose a set of discrete e's. That function f(e) = r{E+ie) can be analytically 
continued to e =0. Pade approximant is well known as a method for the analytical 
continuation. It is useful when coefficients of the Taylor series are given. On the 
contrary, the functional values at several points are given in the present case. We 
employ the point method 8 ) and consider the continued fraction 



1 I Ql(e-ei) 
1 1+- 



_ t(E + eii) ai(e - ei) a 2 (e - eg) , 

1+ 1+ 1+ 1 ' ' 

The coefficients cij are easily deduced from the discrete set of f(e)-values: 

"' '/ '/-I " 1+ ' ' 1+ 

Q l( e '+l- g l) i (2 . 5) 

+ i- [T{E + ie l )/r{E + ie l+1 )Y { ° j 
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e [MeV] 


r{E + ie) [fm 2 ] 


0.1 


-0.11701981630329802 


-0.60857451562152798 i 


0.2 


-0.11924707235802381 


-0.61103993137637069 i 


0.3 


-0.12151390252710244 


-0.61348507952563247 i 


0.4 


-0.1238200220019596 


-0.6159089815519192 i 


0.5 


-0.12616511343686601 


-0.61831066360816722 i 


0.6 


-0.12854882700065837 


-0.62068915788211088 i 


0.7 


-0.13097078049779975 


-0.62304350396337027 i 


0.8 


-0.13343055955913999 


-0.6253727502090104 i 


0.9 


-0.13592771790246552 


-0.62767595510338259 i 


1.0 


-0.13846177766262588 


-0.62995218860804847 i 





-0.1148323866513037 


-0.60608981409107898 i 


exact 


-0.11483238665130231 


-0.60608981409107832 i 



Table I. The propagator r from Eq(2-3) for different e - values. The value for e=0 is evaluated 
by the point method K The exact value is obtained from Eq.(2-3). The parameters are A=- 
0.5592[fm 2 ], / 3=1.13[fm- 1 ] and m _1 =41.47[MeVfiii 2 ]. The energy is 10.0 MeV. 



e [MeV] 


Too [fm 2 ] 


T 02 [fm 2 ] 


T22 [fm 2 ] 


0.4 


0.074878-0.38063 i 


-9.7683 xlO 


~ 4 +8. 8418x10" 


3 i 


4.5582x10" 


a 


2.9901 xlO" 4 i 


0.5 


0.072896-0.38404 i 


-8.0207x10" 


4 +9. 1311x10"' 


1 i 


4.5523x10" 


3 


3.1989xl0~ 4 i 


0.6 


0.070821-0.38742 i 


-6.2122 xlO 


- 4 +9.4171xl0" 


3 i 


4.5460x10" 


3 


3.4052 xl0" 4 i 


0.7 


0.068651-0.39077 i 


-4.3437 xlO 


~ 4 +9.6996xl0" 


3 i 


4.5392x10" 


3 


3.6090 xlO" 4 i 


0.8 


0.066387-0.39409 i 


-2.4160 xlO 


~ 4 +9.9782xl0" 


3 i 


4.5319x10" 


3 


3.8100xl0" 4 i 





0.081860-0.36687 i 


-1.6137 xlO 


" 3 +7.6562 xlO 


- A i 


4.5768 xlO 


-3 


-2.1337 xlO" 4 i 


exact 


0.081860-0.36687 i 


-1.6138 xlO 


" 3 +7.6564 xlO 


" 3 i 


4.5767 xlO 


-3 


-2.1338 xlO" 4 i 



Table II. The T-matrix for the 3 Si- 3 Di state and different e-values. The value for e=0 is calculated 
by the point method and compared to the exact number resulting directly from Eq.(2-7). 



and 

ai = [r(E + ^)ME + i e 2 )]-^ ^ 
£2 — ei 

An example is displayed in Table I, which shows the smooth behavior of f (e) sampled 
for discrete e- values according to Eq(2-3). The analytically continued value for e =0 
as given by (2-4) agrees perfectly well with the exact one. 

Now that approach also works perfectly for a modern high precision NN force 
like the CD-Bonn potential 15 -*. As an example we take the coupled channel case 
3 Si — 3 D\ and solve now the two- body Lippmann Schwinger equation 

f°° o 1 
Tn(p',p) = V n (jp',p) + ]T / dp"p" 2 V n , l (jp'y) tjtTMp",p) (2-7) 

numerically for a set of e- values. We use the same method. The continued fraction 
is now in the matrix 

f W . T<p, p; E + K ) = ^£ + ^, («-«,) „ (M) 

where Ai are components of 2x2 matrices. We arrive at the results displayed in Table 
II. There we show the on- shell T- matrix elements Tqq (p, p) , T02 (p, p) , T21 (p, p) for 
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p = y/mE, E =10.0 MeV and different e-values. Again we see a smooth dependence 
on e. The resulting matrix from the point method elements (e=0) are compared to 
the values evaluated by a standard subtraction method in Table II. The agreement 
is again perfect. 

That expansion works equally well for off- shell T- matrices lim^o T(p',p; E+ie) 
where p / p' and E not connected to p or p' . 

§3. An example for three-nucleon scattering 

We now proceed to three-nucleon scattering. Several methods have been de- 
veloped to treat singularities for the Faddeev equations. The contour deformation 
method has been used as one of the methods for avoiding complicated singularities. 
Here we adopt a separable expansion method 19 )' 20 > for calculating the two-body 
t-matrix. The resulting three-body equations become coupled integral equations in 
one variable and the contour deformation can be used as a tool for avoiding singular- 
ities. We compare the three-nucleon scattering amplitudes obtained by the present 
method with those from the traditional contour deformation method. 

The three-body equations we solve are of the following form 

X ( c P( q , q ';E) = Z { c P( q , q ';E) 

+ E rq' ,2 d q / 'Z^\q,q , '-,E)T(E- q ^/2^fJ( q '',q'-,E) (3-1) 
~, Jo 

Eqs. (3T) are a generalization of the Amado-Mitra-Faddeev-Lovelace(AMFL) equa- 
tions , which are three-body equations based on a simple rank-1 separable potential. 
The AMFL equations are described in standard textbooks for the quantum mechan- 
ical three-body problem and the generalization thereof in ref 21 - 1 . 

The three-body amplitude describing elastic scattering is the on-shell n-d am- 
plitude Xc iCj (<?0) Qo', E) with E = %/2/j, and n-d channels q, Cj allowed for the total 
spin and parity J n . The number of n-d channels are two for J 71 " = (l/2) + , (1/2) ~ or 
three for others. The contour deformation utilizes the analytical continuation in the 
momenta q,q',q", while the present method is based on the analytical continuation 
in the energy E. In both cases, the various singularities in Eq.(3-1) are avoided. 

We solve equ(3.1) with the present method at 10 MeV laboratory energy. At 
this energy, it is known that the breakup cross section is already about 10-20 % of 
the total cross section. Thus the moving singularity is quite important. 

In Table III, we show a result for the three-nucleon scattering amplitude for 
J w = (l/2) + based on the AV14 potential 22 ) . Two channels, denoted by d and q 
and reflecting spin states, are coupled for the n-d elastic scattering amplitude. The 
energy is chosen as E = Eq + ie, where Eq = (/q/2^ and go is the initial momentum. 
In contrast to the simpler two-body singularity treated in section 2 , where real 
epsilons were chosen, we use now also complex epsilons. This allows to have more 
complex energy points in the neighborhood of Eq. Also we kept the imaginary parts 
of the energy always larger than 1 MeV ( this is in relation to the chosen Eq=W 
MeV) which guaranteees a relatively smooth integrand and avoids problems in the 
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ie [MeV] 


M 5 ] 


x ( P [ fm5 ] 


Xip [fm 5 ] 


0+1.0 i 
1.0+1.5 i 
0.0+2.0 i 
0.5+1.5 i 
-0.5+1.5 i 
0.5+1.0 i 
-0.5+1.0 i 
0.5+2.0 i 


0.5634 -1.7745 i 
0.6669 -1.8203 i 
0.7608 -1.8267 i 
0.6513 -1.9441 i 
0.6937 -1.7363 i 
0.4892 -1.9070 i 
0.6265 -1.6940 i 
0.7819 -1.9252 i 


-1.3755 xl0~ 2 +1.4227xl0~ 2 i 
-1.9726 xl0- 2 +8.4048xl0" 3 i 
-2.5101 xl0- 2 +2.3236xl0" 3 i 
-2.5979 xl0~ 2 +1.4182xl0 -2 i 
-1.4067 xl0~ 2 +2.8661xl0~ 3 i 
-1.9834 xl0~ 2 +2.0626xl0~ 2 i 
-8.3854 xl0~ 2 +8.0987xl0 _3 i 
-3.1597 xl0- 2 +7.5166xl0" 3 i 


0.2618+4.5351 xl0~ 2 i 
0.2451+8. 7036xl0~ 2 i 
0.2198+0.1208 i 
0.2878+0.1136 i 
0.2109+6.8861 xl0~ 2 i 
0.3135 +5.9376xl0~ 2 i 
0.2225+ 3.6096xl0~ 2 i 
0.2514+0.1543 i 




C. D. 


0.4248 -1.5440 i 
0.4244 -1.5441 i 


2.5154 xl0 _5 + 2.4663xl0 -2 i 
4.9443 xl0 _5 + 2.4784xl0~ 2 i 


0.2586 -4.5922 xlO" 2 i 
0.2585 -4.5875 xlO -2 i 



Table III. T-matrix for three-body scattering for different complex e- values and J w = (l/2) + . The 
value on the real axis (ie=0) is calculated by the point method and compared to the contour 
deformation number (C. D.) resulting directly from Eq.(3T). The channels d and q correspond 
to spin doublet and quartet, respectively. 



ie [MeV] 


X ( /P [fm 5 ] 


X ( /P [fm 5 ] 


X ( q p [fm 5 ] 


0+1.0 i 
1.0+1.5 i 
0.0+2.0 i 
0.5+1.5 i 
-0.5+1.5 i 
0.5+1.0 i 
-0.5+1.0 i 
0.5+2.0 i 


0.2117 -0.3215 i 
0.3032 -0.2590 i 
0.3575 -0.1809 i 
0.3999 -0.3379 i 
0.2454 -0.1941 i 
0.2813 -0.4530 i 
0.1778 -0.2291 i 
0.4549 -0.2156 i 


-0.2277 -0.2446 i 
-0.1941 -0.2410 i 
-0.1650 -0.2328 i 
-0.1851 -0.2748 i 
-0.1984 -0.2116 i 
-0.2250 -0.2838 i 
-0.2271 -0.2121 i 
-0.1526 -0.2606 i 


-0.6525 -0.6931 i 
-0.5818 -0.6991 i 
-0.5161 -0.6986 i 
-0.5806 -0.7702 i 
-0.5775 -0.6340 i 
-0.6601 -0.7722 i 
-0.6406 -0.6249 i 
-0.5111 -0.7635 i 




C. D. 


-4.5247xl0 -2 -0.3098 i 
-4.5393xl0 -2 -0.3098 i 


-0.3013 -0.2316 i 
-0.3016 -0.2312 i 


-0.7976 -0.6440 i 
-0.7975 -0.6437 i 



Table IV. Same as table III except for J" = (1/2)". 

convergence of the integration. As a consequence of that large imaginary parts of the 
discrete set of energies the amplitudes shown in Table 3 differ considerably from the 
amplitudes, CD., on the real axis as obtained with the contour deformation method. 
Nevertheless the point method applied to those amplitudes leads to good converged 
results which are very close to the amplitudes CD. 

Table IV shows the amplitudes for J 11 = (l/2)~. We see again the same good 
convergence to the results achieved by the contour deformation method. 

§4. Outlook 

We demonstrated the easiness by which the singularities in the Lippmann Schwinger 
equation for NN scattering can be avoided. We used a set of complex energies above 
the positive real axis and based on that performed an analytical continuation to 
the real axis. This analytical continuation was performed by the point method ( a 
continued fraction expansion). 

In case of the three-body Faddeev equations one encounters in the integral kernel 
two types of singularities related to the nucleon-deuteron and three-nucleon cuts. 
The first one is generated by the off- shell NN T- matrix and the second one by the 
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free three-nucleon propagator. As shown in section 3 the complex energy method 
can again be applied. The resulting accuracy is quite satisfactory and the resulting 
amplitudes agree very well with the ones from the contour deformation method. 

This method also allows to search for resonances below the real axis. In a 
forthcoming paper we will display that this method is applicable to even more than 
three particles like the complex four- body break- up process , where no calculations 
have been performed so far. 
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